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Abstract: Based on experiments of the compaction of granular materials under periodic 
shear of a packing of glass beads, a minimal model for the dynamics of the packing density 
as a function of time is proposed. First, a random "energy landscape" is created by a random 
walk (RW) in energy. Second, an ensemble of RWs is performed for various temperatures 
in different temperature-time sequences. We identify the minimum (mean) of the energy 
landscape with the maximum (random) density. The temperature scaled by the step-size 
of the energy landscape determines the dynamics of the system and can be regarded as the 
agitation or shear amplitude. The model reproduces qualitatively both tapping and shear 
experiments. 
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1 Introduction 



The issue of slow dynamics, rare events and anomalous diffusion is object to ongo- 
ing research in statistical physics fl[7|,[l3, 21, 22]. One example for an experimental 
realization of slow dynamics is the compaction of granular materials, being of in- 
terest for both industrial applications and research. Many compaction experiments 
have been carried out in the last decades, see e.g. [fiJHB, 12], but the evolution of 



density with time is far from being well understood. Recent laboratory experiments 
concern pipes filled with granular media and periodically accelerated in order to 
allow for some reorganization [||, |n|,[l8|]. The compaction dynamics was obtained 
to be logarithmically slow and could be reproduced with a simple parking model 
p|,[Tl|]. Alternative experiments concern a sheared block of a granular model ma- 
terial (monodisperse glass spheres) and display a similar dynamics Jig]. Numer- 
ical model approaches, like a frustrated lattice gas (the so-called "Tetris" model) 
[|5, 15, 20] also lead to this slow dynamic behavior, as well as some theory based 



on stochastic dynamics [14]. The fact that a peculiar dynamics is reproduced by so 
many models indicates that it is a basic and essential phenomenon. 

Rather than modeling granular systems in all details, e.g. in the framework of molec- 



ular dynamics or lattice gas simulations [19,20], we will propose a very simple 
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model based on the picture of a random walk in a random energy landscape, for 
a review see [Q], and even simpler than recent, very detailed considerations in the 
same spirit [Q]. Random walks have been examined, for example, on fractals and 
ultrametric spaces, where the continuous time random walk was introduced in order 
to allow a mathematical treatment JlQ|]. If a random walker is situated in a uncor- 
rected, random, fractal energy landscape, the process is called Sinai-diffusion [Qj. 
The aim of this study is to show that the Sinai model is in qualitative agreement with 
the compaction dynamics of granular media. An issue not addressed here is a quan- 
titative adjustment which can be reached, for example, by introducing correlations 
in the EL. 



2 Summary of the experimental results 

The subject of compaction of granular material 12] has been recently revis- 

ited through careful experiments carried out by Knight et al. and Nowak et al. []9l 



17j|18fl. The experiment consists of a vertical cylinder, full of monodisperse beads, 
which is submitted to successive distinct taps of controlled acceleration (vertical vi- 
bration). The measurement of the mean volume fraction after each tap gives a precise 
information about the evolution of the compaction. From the first experiments per- 
formed with taps of constant amplitude, the increase in volume fraction was found to 
be a very slow process well fitted by the inverse of a logarithm of the number of taps. 
Nowak et al. [[n|,[l8|] have then studied the compaction under taps of variable ampli- 
tude, and showed that irreversible processes occur during the compaction. Starting 
with a loose packing, the evolution of the volume fraction is not the same when in- 
creasing the amplitude of vibration as when decreasing. The first branch appears to 
be irreversible, whereas the second is reversible. 

We have recently performed a compaction experiment based on cyclic shear applied 
to an initially loose granular packing A parallelepipedic box full of beads is 
submitted to a horizontal shear through the periodic motion of two parallel walls 
at amplitude # max , see Fig. [j](a). Compaction occurs during this process, leading 
to crystallisation of the beads in the case of a monodisperse material. The control 
parameter in this configuration is the maximum amplitude of shear 6> max (inclination 
angle of the walls). The measurement of the mean volume fraction shows that 
compaction under cyclic shear is a very slow process as in the vertical vibration 
experiments (typically 5 x 10 4 shear cycles or taps). The higher the shear amplitude 
the more efficient is the compaction (shear amplitudes up to # max < 12.5° were 
examined), when starting with a loose packing of the same initial volume fraction. 
More surprising results arise when the packing is submitted to a sudden change in 
shear amplitude. We have observed that a "jump" in volume fraction occurs which 
is opposite and proportional to the change in 6> max . Sudden increase (resp. decrease) 
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Figure 1 : (a) Sketch of the experimental setup, (b) Evolution of the volume fraction 
as function of shear cycles n: a saturated state is seemingly reached after 6 x 10 4 
cycles with an initial shear angle of 5 degrees; a positive jump in volume fraction is 
observed when the angle is decreased to 2.5 degrees and applied for 5 x 10 3 cycles; 
a negative jump is observed when the angle is increased to 10.2 degrees. 

in the shear amplitude decreases (resp. increases) the volume fraction (Fig. []]). The 
response is very rapid (less than 20 cycles) and quasi-independent of the state before 
the angle change. For more detailed experimental results see JI^]. 



3 The Model 

A naive picture that evolves out from the experimental observations is the analogy 
between the packing of beads and a thermal system seeking for a minimum of energy 
in a very complex potential-energy landscape [Q]. Due to some agitation (shear) 
compaction occurs and the total potential energy of the packing decreases. Starting 
from a loose packing at n = 0, the vibrational or shear excitation can be seen as the 
analog of the temperature in the sense that the excitation allows for an exploration 
of the phase space. In a granular packing of monodisperse spheres, the absolute 
minimum of the energy is obtained for a perfect (fee) or (he) crystal with volume 
fraction <p = 0.74. If the energy landscape is complex with a lot of different scale 
valleys or hills, one can understand that an efficient, fast compaction will be obtained 
with high temperatures: the system is then able to escape the deep local valleys 
and find a valley with lower potential energy. A decrease of the temperature is then 
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needed to explore local fine-scale minima. The goal of this paper is simply to explore 
this idea by studying the dynamics of random walkers on a random landscape (the 
Sinai model), and to show that this very naive picture gives results in qualitative 
agreement with the experimental observations. 

Our model is based on the assumption that all possible configurations of a granulate 
in a given geometry can be mapped onto an "energy landscape" (EL). Since a simple 
two-level system (as used to model the dynamics in simple glasses) does not lead to 
the experimentally observed phenomenology, we assume a fractal energy landscape 
created by a random walk in energy with phase space coordinate x. A typical EL 
with energy V(x) is schematically shown in Fig. ^. The stepsize in energy is AV, 
the mean of the energy landscape is V mean and its absolute minimum is V m - m . Here, 
the EL is symmetric to its center, in order to allow for periodic boundary conditions 
in x. Given some EL, the granulate is now modeled as an ensemble of random 
walkers diffusing on the EL, with a temperature Trw. The analog to the density of 
a granular packing is the rescaled energy of an ensemble of random walkers in the 
energy landscape 

E — V m i n 

v = l ~v — ' (1) 

' mean * mm 

which we denote as density in the following. For the (random) initial configuration, 
the energy will be E V mcan so that v = 0; for the close-packing configuration, 
i.e. all RWs are in the absolute minimum, one has E ps V m - m and thus v = 1. The 
energy landscape has the size L, the ensemble of RWs consists of R random walkers 
and S is the number of steps performed. 

Since the energy of the RWs in the EL corresponds to the potential energy of the 
packing, we interpret the (constant) stepsize AV (used here as one possibility for 
the construction of the EL) as a typical activation energy barrier. The maximum of 
V corresponds to a random loose, local packing density, the mean to the (initial) 
random close packing, and the absolute minimum to the hexagonal close packing. 
For one RW, the probability to jump within time At from one site X{ to its neighbor- 
site Xi±\ is 

p±(xi) =mm[l,ex.p(-A ± (xi)/T RW )} , (2) 

with A±(xi) = V{xi±\) — V(xi). Since AV is the only energy scale of the system, 
we define the dimensionless energy steps 5± = A±(xi) / AV and the dimensionless 
temperature T = Trw I A V. Written in dimensionless parameters, the jump prob- 
abilities are thus p l ± = p±(x,f) = min [1, exp(— 6^,/T)] , so that a particle always 
jumps downhill (5± < 0), but jumps uphill only with a probability eo = exp(— 1/T) 
(for 6± > 0), at finite temperature. The limits T — > and T — > oo correspond thus 
to immobile particles or to a homogeneous RW, respectively. The discrete master 
equation for the probability density n l {t), to find a particle at time t at site i of the 
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Figure 2: Schematic plot of a typical energy landscape as used for the simulations. 
The insets show different zoom levels, in order to get an idea of the local situation. 
For an explanation of the symbols, see the text. 



EL, is 



n\t + At)-n\t) 



P++P- 



n\t) + p!j"V +1 (i) + p'+V- 1 ^) (3) 



and is (straightforwardly) simulated with R = 200 random walkers in an energy 
landscape with L = 5000 sites, if not explicitly mentioned. The time interval At 
corresponds to one Monte-Carlo step but not to one shear-cycle n. For the sake of 
simplicity, we measure x in units of the distance between neighboring sites Ax = 
Xi + i — Xi = 1, and time in units of At. The diffusion constant of a homogeneous 
random walk (p± = 1/2 or T — > oo) is thus D = Ax 2 / At = 1. For a constant 
occupation probability (initial density n c = n l (t = 0) = R/L = const.), one can 
extract the diffusion constant D c as function of the temperature, since p± = 1/2 and 
p± = eo/2 occur with equal probability, so that 



DciT) 



1 + e 1 + exp(-l/T) 



(4) 



In Fig. ||, D c = R2(t)/Vi is plotted against T after different times t, with i?2(t) = 
((x(t) — xo) 2 }. For large T, the system does not feel the EL and behaves like an 
ensemble of homogeneous RWs, whereas its behavior becomes subdiffusive after 
several steps for small T when the EL is explored. In a situation with T 0, after 
a rather short transient, all RW will occupy a local minimum, a valley (V), wheras 
the unstable local configurations hill (A), or left- (/) and right-slope (\) cannot be 
occupied, if T is small enough, the random walkers stay trapped. Note that this 
statement is strictly true for T = or when a fixed number of shear-cycles or taps 
is implied. In the Sinai diffusion model, for a finite system, the RW will always find 
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Figure 3: (a) Density plotted against time t (in units of Monte-Carlo steps) for a 
simulation with Tq = 1.1, L = 5000 and R = 200. The temperature-sequence is 
indicated by the dashed line, (b) Diffusion constant D c = R,2(t)/y/i at time t as 
function of temperature T for systems with constant (initial) density. The solid line 
shows Eq. (|j). 



the global minimum - the temperature only determines the time-scale of this process 
Since the duration of an experiment is limited, the global minimum cannot be 
found by all particles if the phase space volume L is large enough. 



4 Results and Discussion 



In parallel to the tapping experiments by Nowak et al. [|18fl, we present simulations 
of our model using a periodic time series for the temperature. The temperature is 
kept constant for S steps and then increased by AT = 0.1, where it is kept again 
constant for S steps. T is initially zero and then raised up to T = 2. From this state, 
T is decreased to zero and the loop is repeated seven times. In Fig. ^, the simulation 
results are displayed for different S as given in the inset. In the initial branch with 
increasing T, the density increases and slowly decreases for large T. This branch is 
irreversible, but the periodic loops show almost reversible behavior. For very short 
loops (small S), the density continuously increases, for longer loops the behavior 
of the system is reversible. Note that the system also shows hysteretic behavior, the 
density at decreasing T is below the density at increasing T. 

In summary, we presented a very simple model for the dynamics of the compaction 
of granular media due to an external agitation. Our model is not as detailed as others 
(parking lot or frustrated lattice-gas models), but it is extremely simple and still 
shows qualitative agreement with two different types of experiment. Its simplicity 
allows for future analytical treatments. However, such a simple model arises, besides 
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Figure 4: Rescaled density v plotted as a function of the temperature T for a periodic 
temperature sequence. The arrows indicate the direction of the loops. 



many others, two major questions: (i) Is it possible to link the real configuration 
phase space with the an energy landscape? (ii) Is it correct to do the analogy between 
the experimental external excitations like shear or tapping with a temperature? 
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